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Deciphering and understanding the small-scale magnetic activity 
of the quiet solar photosphere should help to solve many of the key prob- 
lems of solar and stellar physics, such as the magnetic coupling to the 
outer atmosphere and the coronal heating 1 ' 2 ' 3 . At present, we can see 
only ~1% of the complex magnetism of the quiet Sun 1 ' 4 ' 5 ' 6 ' 7 , which high- 
lights the need to develop a reliable way to investigate the remaining 99%. 
Here we report three-dimensional radiative tranfer modelling of scatter- 
ing polarization in atomic and molecular lines that indicates the presence 
of hidden, mixed-polarity fields on subresolution scales. Combining this 
modelling with recent observational data 8 ' V } 11 we fi nc i a ubiquitous tan- 
gled magnetic field with an average strength of ~130 G, which is much 
stronger in the intergranular regions of solar surface convection than in 
the granular regions. So the average magnetic energy density in the quiet 
solar photosphere is at least two orders of magnitude greater than that 
derived from simplistic one-dimensional investigations 12 ' 13 , and sufficient 
to balance radiative energy losses from the solar chromosphere. 

Most of our present empirical knowledge of solar surface magnetism stems 
from the analysis of the light polarization that the Zeeman effect induces in spectral 
lines 1 ' 4 ' 5 ' 6 ' 7 . For example, the circular polarization signals of which synoptic magne- 
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tograms are made reveal the existence of an irregular network of spatially unresolved, 
intermittent flux patches of kilogauss field concentrations outlining the boundaries of 
the giant velocity cells of the supergranulation 1 . Such routine magnetograms give the 
(wrong) impression that the cell interiors (referred to here as internetwork regions) 
are non-magnetic, and it is important to note that such regions cover most of the 
solar 'surface' at any given time during the solar magnetic activity cycle. However, 
high-spatial resolution magnetograms show a multitude of mixed magnetic polarities 
within the internetwork regions, with most of the detected flux located in the inter- 
granular lanes of the solar granulation pattern where the plasma is downflowing 4,5 ' 6,7 . 
The mean unsigned flux density turns out to be ~10 G when choosing the best com- 
promise between polarimetric sensitivity and spatio-temporal resolution that is at 
present feasible. When such polarization signals are interpreted taking into account 
the fact that the magnetic field is not being spatially resolved (for example, by as- 
suming that one or more magnetic components are coexisting with a 'non-magnetic' 
component within the spatio-temporal resolution element of the observation), it is 
then found that the filling factor of the magnetic component(s) is ~1%. The 'problem' 
with the Zeeman effect is that the amplitudes of the measured polarization signals 
are the smaller the greater the degree of cancellation of mixed magnetic polarities 
within the spatio-temporal resolution element of the observation. Therefore, lack of 
detection does not necessarily imply absence of magnetic fields. 
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Fortunately, scattering processes in spectral lines produce linear polarization 
signals, whose amplitudes are efficiently modified in the presence of tangled magnetic 
fields of strength 5 H ~l-137x 10 _7 /(tii fc g L ) (with B H expressed in gauss, and where 
iiif e and g L are respectively the lifetime in seconds and Lande factor of the upper level 
of the transition). This so-called Hanle effect 14 ' 15 has the required diagnostic poten- 
tial for investigating spatially unresolved, 'hidden' mixed-polarity magnetic fields in 
the solar atmosphere 16 . The problem is how to apply it to obtain reliable information 
given that Hanle-effect diagnostics of such 'turbulent' fields relies on a comparison 
between the observed scattering polarization and that corresponding to the zero-field 
reference case. It has been pointed out correctly 12 ' 13 that a suitable spectral line for 
Hanle-effect diagnostics of 'turbulent' photospheric fields is that of Sr I at 4607 A. 
However, the simplified approach of assuming that the highly inhomogeneous and 
dynamic solar photosphere can be well represented by a one-dimensional (ID) and 
static model is very unreliable, because of the need to use the free parameters of 'clas- 
sical' stellar spectroscopy (that is, micro- and macroturbulence for line broadening), 
which have a serious impact on the calculated polarization amplitudes 17 . We have 
shown 17 that the particular ID approach that has been applied 12,13 yields artificially 
low values for the strength of the 'turbulent' field -that is, between 20 and 10 gauss 
as shown by the dashed lines of Fig. 3 in Ref. 17. 
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In order to improve the reliability of diagnostic tools based on the Hanle 
effect, we have developed a novel approach based on multilevel scattering polarization 
calculations in three-dimensional (3D) models of the solar photosphere, which we have 
obtained from realistic hydrodynamical simulations of solar surface convection 18 . This 
has allowed us to obtain the linear polarization amplitudes that scattering processes 
in the inhomogeneous solar photosphere would produce in the Sr I 4607 A line if there 
were no magnetic field. We point out that our synthetic intensity profiles (which take 
fully into account the Doppler shifts of the convective flow velocities in the 3D model) 
are automatically in excellent agreement with the observations when the meteoritic 
strontium abundance is chosen. However, we find that the calculated fractional linear 
polarization is substantially larger than the observed one, thus indicating the need 
to invoke magnetic depolarization. 

As seen in Fig. 1, our radiative transfer simulations of the Hanle effect 
in the Sr I line show that a volume-filling, microturbulent magnetic field of about 
60 gauss leads to a notable agreement with the observed Q/I (where Q and / are 
Stokes parameters). In Fig. 1 there is a clear indication that the strength of the 
'turbulent' field required to explain the Q/I observations decreases with height in the 
atmosphere, from the 70 gauss needed to explain the observations at fi = 0.6 to the 
50 gauss required to fit the observations at fj, — 0.1. This corresponds approximately 
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to a height range between 200 and 400 km above the solar visible 'surface'. 

If the hidden magnetic field of the quiet solar photosphere really had a 
strength of 60 G at all points in the solar photosphere, then the mean magnetic 
energy density (E m =< B 2 > /8n) would be E m ttl40 erg cm -3 , which is an or- 
der of magnitude smaller than that corresponding to the kG fields of the network 
patches. However, both magnetohydrodynamical simulations 19 ' 20 and observations 
of the Zeeman effect 4,21 ' 7 indicate that the photospheric plasma of the quiet Sun has 
a continuous distribution of field strengths, in such a way that the weaker the field 
the larger the probability of finding a magnetic strength between B and B + dB when 
no distinction is made between granular and intergranular regions. As shown by the 
black dashed-dotted line of Fig. 1, if we assume that the probability distribution 
function (PDF) has an exponential shape (e~ B / B " /B ), we then find that B 130 
gauss yields a fairly good fit to the observed fractional linear polarization. In this 
much more realistic case £ m ^1300ergcm~ 3 (that is, V< B 2 > xs 180 G), which is 
about 20% of the averaged kinetic energy density produced by convective motions 
at a height of 200 km in the 3D photospheric model. This result and the fact that 
the observed Q/I does not seem to be modulated by the solar cycle (see Fig. 1) 
suggests that a small-scale dynamo associated with turbulent motions within a given 
convective domain of ionized gas 19 ' 20 plays a significant role for the quiet Sun mag- 
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netism. The total magnetic energy stored in the internetwork regions is now slightly 
larger than that corresponding to the kG fields of the network patches. It is also 
of interest to point out that our empirical PDF with < B > xs 130 G produces a 
negligible Zeeman broadening of the emergent spectral line profiles, which is fully 
compatible with the reported empirical constraints 1 . However, in contrast to what 
might be expected from the conclusions of a recent paper 22 , the stretched exponential 
PDF that characterizes the surface distribution of magnetic fields 'predicted' by the 
most recent numerical experiments of turbulent dynamos 19 , which correspond to a 
magnetic Reynolds number of about 1000, implies a Hanle-effect reduction of the 
scattering polarization amplitude of the Sr I 4607 A line that is significantly smaller 
than what is needed to explain the observations. 

We do not have a sufficient number of observational constraints to be able to 
determine the exact shape or the detailed spatial variation of the solar PDF. However, 
we can conclude that the 'turbulent' field is organized at the spatial scales of the solar 
granulation pattern, with relatively weak fields above the granules and with much 
stronger fields above the intergranular lanes (see Fig. 2). As shown in Fig. 2 legend 
most of the volume of the upflowing granules where the observed C 2 line polarization 
originates is occupied by magnetic fields much weaker than what is needed to explain 
the observed depolarization in the Sr I 4607 A line. Therefore, most (but not all) 
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of the observed strontium line depolarization must be produced by relatively strong 
and tangled fields in the intergranular regions. 

The contribution of the intergranular plasma to the observed scattering po- 
larization in the strontium line is very likely to be close to the regime of Hanle 
saturation, which for this spectral line occurs for magnetic strengths B>300 G. This 
helps us to understand why with low spatial and temporal resolution it is not easy 
(but definitely not impossible) to detect spatial fluctuations in the observed Q/I 
(Fig. 1). We have carried out calculations with several plausible PDFs, assum- 
ing in all cases that the angular distribution of the field vectors is isotropic and 
microturbulent for every single field strength. For instance, we have assumed an 
exponential PDF (PDF G = e~ B/B z/B g ) for the upflowing regions with a free pa- 
rameter B g to be determined from the observations and a given maxwellian function 
(PDFi = 2.38 x 1CT 8 -B 2 exp[— B/456} 2 ) for the downflowing plasma obtained from a 
best fit to the strong field part of the intergranular histogram of the Zeeman splittings 
observed in near-IR lines of neutral iron 7 . Note that this maxwellian PDF implies 
that the filling factor of kG fields is ~2% -that is, it does not exclude the possibility 6 
of small-scale kG fields in the internetwork regions. With this 'educated' choice of 
plausible PDFs the best agreement between theory and observations is obtained for 
Bg&15 G, which is consistent with the results of our investigation of the Hanle ef- 
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feet in C2 lines. We find that most of the magnetic energy turns out to be due to 
rather chaotic fields in the intergranular plasma with strengths between the equipar- 
tition field values and ~1 kG. The total magnetic energy stored in the internetwork 
regions is now substantially larger than that corresponding to the kG fields of the 
supergranulation network. 

Our conclusion that most of the volume of the granular regions is occupied 
by very weak fields seems to be compatible with the constraints imposed by our 
present understanding 23 of the 'enigmatic' scattering polarization observed in the Na 
I D x line, whose line-core radiation originates a few hundred kilometers higher in the 
solar atmosphere. As demonstrated in Ref. 23, the linear polarization signature of 
the scattered light in the Di line is only completely suppressed in the presence of 
magnetic fields larger than 10 G, regardless of their inclination. Note that this result 
is significantly different to what had been concluded in an earlier paper 24 . In any 
case, it is important to point out that the observed linear polarization in the sodium 
Di line still remains enigmatic because nobody has yet been able to model both the 
amplitude and shape of the observed Q/I profile 25 . 

Our empirical findings may have far-reaching implications in solar and stellar 
physics. The hot outer regions of the solar atmosphere (chromosphere and corona) 
radiate and expand, which takes energy. By far the largest energy losses stem from 
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chromospheric radiation with a total energy flux of 26 ~10 7 erg cm~ 2 s -1 . Con- 
sidering our most conservative estimate for the magnetic energy density -that is, 
£' m ~140ergcm~ 3 - we obtain an energy flux similar to the above-mentioned chro- 
mospheric energy losses when using either the typical value of ~lkms~ 1 for the 
convective velocities or the Alfven speed (va = B /(4-rcp) 1 / 2 , with p the gas density). 
In reality, as pointed out above, the true magnetic energy density that at any given 
time during the solar cycle is stored in the quiet solar photosphere is very much larger 
than 140 erg cm" 3 . Only a relatively small fraction would thus suffice to balance the 
energy losses of the solar outer atmosphere. 

Equally interesting is our result that < B > 130 G, which indicates that 
the unsigned magnetic flux density in the quiet solar photosphere is substantially 
large. In fact, it has been suggested recently that the dynamic geometry of the 
magnetic connection between the photosphere and the corona may sensitively depend 
on the amount of magnetic flux that exists in the internetwork regions 3 . The hidden 
flux that we have diagnosed turns out to be organized at the spatial scales of the 
solar granulation pattern, with much stronger fields in the intergranular regions. 
Therefore, if the magnetic field of the 'quiet' solar chromosphere is found to be 
more complex than simply canopy-like, because of the influence of the intergranular 
magnetic fields which are continuously moving around on temporal scales of minutes, 
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then the possibilities for magnetic reconnection 27 and energy dissipation in the solar 
outer atmosphere would be much enhanced. The 'hidden' magnetic fields that we 
have reported here could thus provide the clue to understanding how the 10 6 K solar 
corona is heated. 
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Figure 1. Spectropolarimetric observations versus 3D modelling of 
the Hanle effect. 

This figure shows the center-to-limb variation of the fractional linear po- 
larization at the core of the Sr I 4607 A line after substraction of the continuum 
polarization level. (The Stokes / and Q parameters are defined in Ref. [1]; jj, = cos 9, 
with 9 the angle between the solar radius vector through the observed point and the 
line of sight). Open circles, various observations 8 taken during a minimum period of 
the solar cycle (in particular, during September-October 1995). The remaining sym- 
bols correspond to observations taken during the most recent maximum period of the 
solar cycle. Diamonds, observations 9 obtained during May 2000; 'crosses' and 'plus' 
symbols, observations taken during August 2000 (see Ref. 10) and December 2002 (V. 
Bommier 2002; personal communication), respectively; filled circles, observations ob- 
tained during September 2003 at the Istituto Ricerche Solari Locarno (Switzerland) 
in collaboration with M. Bianda. Coloured lines, the results of our 3D scattering 
polarization calculations in the presence of a volume-filling and single-valued micro- 
turbulent field with an isotropic distribution of directions below the mean free path 
of the line photons (from top to bottom: 0, 5, 10, 15, 20, 30, 40, 50, 60, 80, 100, 150, 
200, 250 and 300 gauss). We have solved the relevant equations 28 via the application 
of efficient radiative transfer methods 29 , and using realistic collisional depolarizing 
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rate values 12 , which turn out to be the largest rates among those found in the liter- 
ature. Note that there is no evidence of a serious modulation of the strength of the 
'turbulent' field with the last solar activity cycle, and that the best average fit to the 
observations is obtained for 60 G. The black dashed-dotted line indicates the result- 
ing Q/I line-core amplitudes for the case of a single exponential PDF (e~ B ^ B ° /B ) 
with B = 130 G, which implies < B > = J BPDF(B) dB = 130 G. 
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Figure 2. Evidence for a 'turbulent' field organized at the spatial 
scales of the solar granulation pattern. 

This figure demonstrates that the scattering polarization observed in molec- 
ular lines 11 comes mainly from the upflowing regions of the 'quiet' solar photosphere. 
The solid line contours delineate such upflowing regions at two heights in the 3D 
photospheric model, which indicate the approximate atmospheric region where the 
observed C 2 line radiation originates. The two top panels show the horizontal fluctu- 
ation of the 'degree of anisotropy' (A = Jq/ Jq) 28,15 of the solar continuum radiation 
at 5000 A . (Note that there is a very strong correlation between the upflowing regions 
and A). The two bottom panels visualize the corresponding horizontal fluctuation in 
the number density (A/") of C 2 molecules. (Note that there is a significant correlation 
between the upflowing regions and Af). Since in the solar atmosphere the scattering 
polarization in weak molecular lines is proportional to both A and Af/r) c (where r\ c is 
the background continuum opacity, which is anticorrelated with Af in such range of 
heights) we conclude that only the upflowing regions make a significant contribution 
to the observed molecular scattering polarization. Therefore, we can obtain informa- 
tion on the distribution of magnetic fields in such (granular) upflowing regions of the 
solar photosphere via theoretical modelling of the observed scattering polarization in 
suitably chosen pairs of C 2 lines 30 . We have calculated the magnetic sensitivity of the 
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observed C2 lines by taking into account all the relevant optical pumping mechanisms 
in a very realistic multilevel model for C 2 . Thanks to this Hanle-effect investigation, 
we have been able to conclude that in the granular regions the 'turbulent' field is 
much weaker than what is needed to explain the observed depolarization in the Sr I 
4607 A line. For this reason, the intergranular regions must be pervaded by relatively 
strong tangled fields capable of producing most of the observed depolarization in the 
strontium line. 
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